AM, 
E D 学 报 http : // www. insect. org. cn 











基于 mtCOI X Al Fr 71] ^J rh Es] 2 EXER ig 
Ae EE 贵 传 多 样 性 分 析 





An 


HART, REOR EUR, EEE, PERN, A nen 


(1. 长 江 大 学 农学 院 , 长 江 大 学 昆虫 研究 所 , WIIN 434025; 2. 中 国 农业 科学 院 植物 保护 研究 所 ， 
植物 病虫害 生物 学 国家 重点 实验 室 , 北京 100193 ) 






































WE: [An] a Liriomyza chinensis 是 新 蒜 类 蔬菜 上 的 重要 经 济 害 虫 , 在 我 国 广泛 分 布 。 本 
旨 在 冰 析 中 国 茜 斑 潜 晶 的 地 理 种 群 遗 传 分 化 。【 方 法 】 以 中 国 8 省 12 个 不 同 地 理 种 群 的 253 
头 芋 斑 潜 蝇 为 样本 ,测定 其 mtCOI 基因 序列 ;依据 获得 的 mtCOI 基因 序列 ,利用 MEGA7. 0, DnaSP 
6.1 和 Arlequin 3.5 等 软件 对 茜 斑 潜 晶 地 理 种 群 的 遗传 多 样 性 基因 流水 平和 遗传 变异 等 进行 分 
析 。【 结 果 】 在 253 头 个 体 的 759 bp mtCOI 基因 片段 中 ,获得 13 个 单 倍 型 ,各 单 们 型 间 K2P 遗传 距 
离 均 小 于 0.02; 其 中 单 倍 型 Hapl 为 12 个 地 理 种 群 所 共享 ,总 发 生 频率 高 达 81.82% 。 总 群体 单 倍 
型 多 样 性 较 低 ( Hd 20.327) , 核 普 酸 多 样 性 (Pi) A 0.00159 , J£ 3F 88 2E ŽAK) A 1.21011; X 
斑 潜 蝇 总 群体 遗传 分 化 程度 中 等 (ss =0.06971) ,基因 交流 较 充 分 (Nm =3.33629 ) 。 分 子 方差 分 
析 (AMOVA ) 结果 表明 ,遗传 变异 来 源 为 群体 内 ;总 群体 Tajima”s D 检验 值 为 显著 负 值 ;Mantel 检 
测 结 果 说 明 种 群 间 的 遗传 距离 与 地 理 距 离 没 有 相关 性 。【 结论 中 国 芋 斑 漆 蝇 地 理 种 群 的 遗传 多 
样 性 较 低 ,基因 交流 较 充 分 ,遗传 分 化 程度 中 等 , 且 地 理 距离 并 不 影响 其 遗传 分 化 程度 ; 草 斑 潜 晶 总 
群体 在 较 近 的 历史 时 期 未 经 历 明 显 的 种 群 扩 张 和 种 群 增长 。 
关键 词 : 黄斑 潜 蝇 ; 线粒体 COI 基因; 地 理 种 群 ; 遗传 多 样 性 ; 系统 进化 分 析 
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Analysis of the genetic diversity of geographical populations of Liriomyza 
chinensis ( Diptera: Agromyzidae ) in China based on mtCOI gene 


sequence 

ZHONG Yu-Jun'^, DU Su-Jie’, PAN Li-Ting , WANG Yu-Sheng’, WANG Fu-Lian"*, LIU Wan- 
Xue^" (1. Institute of Entomological Science, College of Agriculture, Yangtze University, Jingzhou, 
Hubei 434025, China; 2. State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute 
of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China) 

Abstract: [ Aim] Liriomyza chinensis, a serious economic pest of onions and garlics, is widely 
distributed in China. This study aims to analyze the genetic differentiation among the geographical 
populations of L. chinensis in China. [Methods] A total of 253 individuals from 12 different geographical 
populations from 8 provinces were used as samples, and their mtDNA COI gene was sequenced. The 
genetic diversity, gene flow levels and genetic variation of the geographical populations of L. chinensis 
were analyzed by MEGAT. 0, DnaSP 6. 1 and Arlequin 3. 5 based on the obtained mtCOI gene 
sequences. [Results] From the 253 individuals 13 haplotypes of the mtCOI gene fragment of 759 bp in 
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length were defined, and the K2P genetic distance between haplotypes was less than 0. 02. Haplotype 
Hapl was shared by the 12 geographical populations, and the total occurrence frequency was up to 
81. 8296. The haplotype diversity ( Hd) of all geographical populations was low (0. 327) , the nucleotide 
diversity ( Pi) was 0. 00159, and the average number of nucleotide difference (K) was 1. 21011. The 
total population had generated moderate genetic differentiation (Fs — 0. 06971) , and the gene flow was 
relatively sufficient (Nm 23.33629). The results of molecular variance analysis ( AMOVA) showed that 
the source of genetic variation was within population, and the Tajima’ s D test value was significantly 
negative. Mantel test results showed that the genetic distance was not related to the geographical 
distance. [Conclusion] The geographical populations of L. chinensis in China show low genetic diversity , 
relatively sufficient gene flow and moderate genetic differentiation, and the geographical distance does not 
affect the degree of genetic differentiation between geographical populations. The total population of L. 
chinensis has not experienced significant population expansion and population growth in the recent 
historical period. 


Key words: Liriomyza chinensis; mtDNA COI gene; geographical population; genetic diversity; 


phylogenetic analysis 





ZU XXE 1E B Liriomyza chinensis 3k Jg X00 38 H 
( Diptera) 潜 晶 科 ( Agromyzidae ) 斑 潜 量 属 Liriomyza , 
J&—Rh4E AEE f dE E MB IS SE. d, E EDAUSECKU 
Allium fistulosum 韭菜 4. tuberosum JEX, A. cepa. 等 

















^E, H 2014 年 农业 农村 部 种 植 业 管理 司 的 统计 结 
果 , 我 国 大 萄 种 植 面积 已 逾 57.3 万 hm , 主 产 区 位 
于 又 岭 淮 河 以 北 的 广大 地 区 ,涵盖 华北 平原 、 东 北平 
原 和 黄土 高 原 三 大 主 产 区 ,其 中 华北 平原 的 山东 、 河 

















百合 科 ( Liliaceae ) RJE Allium 植物 ,以 大 营 受 害 最 
(Spencer, 1973; 康乐 ,1996 ) 。 其 幼虫 、 成 虫 均 
可 为 害 , 主 要 以 幼虫 潜 食 叶肉 为 主 ,叶肉 被 取 食 后 ， 
仅 剩 两 层 表 皮 ,呈现 白色 线 状 条 斑 , 形 成 潜 道 , 严 
影响 植物 正常 的 光合 作用 ;并 且 上 肉 成 虫 的 产 卵 器 可 
刺 破 叶片 ,由 此 造成 的 机 械 损伤 为 病原 物 侵 入 提供 
了 重要 途径 ( 潘 秀 美和 姜 官 恒 , 2001; 魏 国 先 等 ， 
2001)。 葱 斑 潜 蝇 广 泛 分 布 于 东南 亚 和 太平 洋 地 
区 ,在 中 国 、 日 本 .韩国 印度尼西亚、 新 加 坡 、 朝 鲜 、 
越南 等 国家 危害 严重 (Tran and Takagi, 2005; Tran 
et al., 2007) ,近年 来 该 虫 也 在 欧洲 部 分 地 区 如 法 
国 \ 德 国 等 地 暴发 成 灾 ( Martynov et al., 2016) 。 该 
虫 在 我 国 大 部 分 地 区 均 具有 和 较 高 的 适应 性 ( 梁 超 峰 
等 , 2018) ,主要 分 布 于 黑龙 江 、 内 蒙古 .辽宁 、 沈 阳 、 
山东 、 河 北 、 河 南 、 山 西 .福建 ,台湾 等 省 区 ( 王 永 卫 





























南 和 河北 三 省 为 我 国 大 萄 的 第 一 大 产 区 ,总 栽培 面 
积 占 全 国 的 36% ,总 产量 占 全国 的 51% (彭帅 等 ， 
2017) 。 山 东 省 作为 我 国营 和 蒜 类 蔬菜 种 植 面积 最 大 
的 省 份 , 2001 年 新 蒜 类 蔬菜 种 植 面 积 约 26.7 万 
hm ,其 80% 的 面积 遭受 莹 斑 潜 蝇 的 危害 ,被 害 株 率 
常 达 40% 以 上 ( 潘 秀 美和 姜 官 恒 , 2001) 。 傅 苏 友 
和 苏 保 乐 (2011 ) 在 调查 山东 潍坊 市 的 葱 斑 潜 晶 发 
生 情 况 时 ,发 现 其 在 大 萄 上 为 害 严重 ,一 株 萄 上 幼虫 
平均 为 7 头 ,被 害 株 高 达 70% ~ 8096 ,严重 的 达到 
100% 。 

目前 有 关 瓯 斑 潜 蝇 的 研究 多 局 限于 生物 学 特 
性 ,发生 规律 等 基础 性 研究 ,而 有 关 葱 斑 潜 蝇 的 种 群 
遗传 多 样 性 、 种 群 扩张 等 方面 研究 还 暂 未 见报 道 。 
新 斑 潜 蝇 作 为 一 种 本 地 害虫 ,相对 于 斑 潜 蝇 属 的 其 
fli A UR, ARE OS PEINE L. huidobrensis 、 美 洲 斑 潜 





















































等 , 1992; 王 进 忠 等 , 2000; Shiao, 2004) , ZR PEX 
晶 的 温度 适应 范围 广 , 单 肉 产 卵 量 大 ,在 我 国 各 地 一 
ERE 4 ~7 代 不 等 ,如 东北 4~5 代 , 山 东 6 ~7 代 ， 
Het dcs AE (王立 霞 等 , 2000; 潘 秀 美和 姜 官 
恒 , 2001) ,有 春季 和 秋季 两 个 发 生 高 峰 期 ,高 峰 期 
易 在 大 萄 种 植 区 肾 发 成 灾 ( 潘 秀美 和 姜 官 恒 , 2000; 
王立 霞 等 , 2000; Tran et al., 2007) 。 

大 葡 作 为 萄 斑 潜 蝇 的 主要 寄主 植物 ,其 大 规模 
化 面积 种 植 一 定 程度 上 有 利于 萄 斑 潜 蝇 的 严重 发 


























ha L. sativae =U E JETER WE L. trifolii 等 ,其 长 期 累 
积 的 遗传 变异 是 否 与 其 严重 危害 存在 联系 ,还 尚 不 
清楚 。 分 子 标记 技术 的 快速 发 展 ,很 大 程度 加 快 了 
内 外 学 者 对 斑 潜 蝇 属 害虫 遗传 分 化 的 研究 。 
Morgan 等 (2000) 首先 对 美国 加 利 福 尼 亚 中 部 和 南 
部 的 南美 斑 潜 蝇 种 群 间 进行 RAPD-PCR 检测 , 发现 
南美 斑 潜 蝇 种 群 间 的 遗传 变异 主要 来 源 于 地 理 分 
布 。 王 莉 萍 等 (2007) 利 用 rzDNA-ITS1 标记 对 我 国 8 
个 不 同 地 理 的 美洲 斑 潜 蝇 种 群 进行 分 析 , 得 出 各 地 
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理 种 群 之 间 的 遗传 分 化 很 小 , 且 遗 传 分 化 很 有 可 能 
与 地 理 分 布 有 关 。 本 研究 选用 线粒体 COI 基因 作 
为 分 子 标记 ,以 我 国 8 省 12 个 不 同 地 区 的 萄 斑 潜 昌 
种 群 为 研究 对 象 ,对 不 同 地 理 种 群 的 遗传 多 样 性 、 遗 
传 分 化 程度 .基因 交流 等 进行 了 分 析 , 以 揭示 其 种 群 
系统 发 育 地 理 格 局 与 演变 机 制 。 人 研究 结果 对 萄 斑 潜 
晶 的 种 群 遗 传动 态 分 析 、 种 群 扩张 机 制 探讨 具有 一 
定 的 指导 意义 ,同时 也 为 制定 葱 斑 潜 蝇 的 防治 策略 
提供 理论 依据 。 
































1 材料 与 方法 


L1 idi 
MERE 2016 - 2018 年 采集 自 中 国 8 
省 12 个 地 点 ( 表 1) ,采用 五 点 取样 法 ,采集 露天 种 
植 的 优势 寄主 植物 大 葱 A. fistulosum BUNT A. 
ascalonicum ,采集 虫 态 为 幼虫 或 师 , 带 回 实验 室 人 饲养 。 
待 虫 源 羽化 后 , 温 于 无 水 乙醇 中 , -40%C 保 存 备用 。 











表 1 葱 斑 潜 蝇 不 同 地 理 种 群 的 采集 信息 


Table 1 Collecting data of different geographical populations of Liriomyza chinensis 



































: 采集 地 点 " 采集 时 间 ( 年 .月 ) 样本 数量 
种 群 经 纬度 寄主 植物 
Sampling Sampling date Number of 
Population Latitude and longitude Host plants 
locality ( year. month ) samples 
吉林 通化 
JLTH 41?46'N, 125?54'E 2018.07 小 萄 Allium ascalonicum 20 
Tonghua, Jilin 
北京 海淀 . 
BJHD 40*01'N, 116?16'E 2017.07/2018. 09 小 萄 A. ascalonicum 21 
Haidian, Beijing 
ans 
NMWH TT 39*48'N, 106^46'E 2017.08 INŽ A. ascalonicum 10 
Wuhai, Inner Mongolia 
河北 保定 
HBBD 38*51'N, 115?35'E 2017.05/2017.07 KAA. fistulosum 20 
Baoding, Hebei 
河北 石家庄 
HBSJZ 37*58'N, 114?26'E 2017.05/2018. 09 /NÆ A. ascalonicum 20 
Shijiazhuang, Hebei 
东 潍坊 
处; 
SDWF 37?49'N, 117?29'E 2017.09 KAA. fistulosum 21 
Weifang, Shandong 
OEN 
AR 烟 r1 N 
SDYT 37°38' N, 120°28' E 2016. 07/2017. 05 KAA. fistulosum 26 
Yantai, Shandong 
ty 
东 济南 
SDJN 37?08'N, 117?06'E 2017.05/2017.07 KARL A. fistulosum 20 
Jinan, Shandong 
东 临 沂 
SDLY 36?18'N, 119*50'E 2018.09 JNA A. ascalonicum 20 
Linyi, Shandong 
甘肃 白银 
GSBY á 36?22'N, 104?18'E 2018.08 KAA. fistulosum 23 
Baiyin, Gansu 
河南 许昌 
HNXC 33*52'N, 113?28'E 2018.01 KAA. fistulosum 24 
Xuchang, Henan 
海南 陵 水 
HNLS 18?28'N, 110?04'E 2018.02 小 萄 A. ascalonicum 28 





Lingshui, Hainan 


1.2 基因 组 DNA 提取 

参照 De Barro 和 Driver( 1997 ) 的 快速 研磨 法 ， 
并 稍 作 改 动 。 具 体 步 又 如 下 :从 标本 管 中 取 单 头 萄 
斑 潜 蝇 雌 成 虫 , 置 于 干净 滤纸 上 , 待 无 水 乙醇 挥发 
后 ,吸取 10 pL 裂解 液 (0. 05 mol/L Tris-HCl, pH 
8.4, 50 mmol/L KCl, 0. 4596 Tween-20, 0.45% NP 
40, 0.296 Gelatin, 80 pg/mL 和 蛋白酶 K) 滴 于 封口 膜 











上 ,用 灭 菌 过 的 0.2 mL PCR 管 管 底 充分 研磨 虫 体 ， 
将 提取 液 吸入 0.2 mL PCR 管 中 。 吸 取 20 uL 裂解 
液 清洗 PCR 管 底 和 封口 膜 2 次 ,最 后 将 清洗 液 吸入 
0.2 mL PCR 管 中 。 将 提取 液 进 行 金 属 浴 , 程序 : 
65C 下 30 min, 25C F2 min, 96% F 10 min, 灭 活 
和 蛋白酶 K, 最 终 获 得 50 uL ARE DNA ,保存 在 
-20C 冰箱 中 备用 。 
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1.3 mtCOL 基因 序列 的 扩 增 

采用 引物 组 合 C1-J-2183 (5'-CAACATTTATTTT 
GATTTTTTGG-3') fll TL2-N-3014 (5'-TCCATTGCAC 
TAATCTGCCATATTA-3') ( Simon et al., 1994) 扩 增 
ZU DER mtCOLI 基因 的 目的 片段 。 引 物 由 上 海 生 
工 生 物 工 程 股份 有 限 公司 合成 。 

PCR 反应 体系 (25 pL): 10 xbuffer( 含 Mg… ) 
2.5 uL, dNTPs(2.5 mmol/L)0.4 uL, 正 反 引物 (10 
pmol/L) 各 0.4 uL, Taq 酶 (2.5 U/uL) 0.2 uL, fX 
板 DNA 1 uL, ddH,O 20.1 uL, PCR 反应 程序 : 
95*C 预 变 性 4 min; 95% 变性 1 min, 52. 5C 退火 50 
s, 72C 延伸 90 s, 35 个 循环 ; 继续 72%C 延伸 10 
min。 取 20 uL 扩 增 产物 在 含有 ColdView 开 的 1% 
琼脂 糖 凝 胶 上 以 110 V 电泳 分 离 30 min, 然 后 以 凝 
胶 成 像 系统 分 析 电 泳 结 果 并 拍照 保存 。 
1.4 PCR 产物 的 纯化 及 克隆 测序 

由 于 部 分 虫 源 直 接 测序 时 出 现 重 奏 峰 ,为 了 确 
保 测序 的 准确 性 ,对 存在 重奏 峰 的 样品 进行 了 克隆 
测序 。 在 紫外 灯 下 切 下 含有 目的 片段 的 明亮 条 带 ， 
使 用 AxyPrep DNA 凝 胶 回 收 试剂 盒 将 PCR 产物 回 
收 纯化 。 将 纯化 产物 与 克隆 质粒 pEASY-T3 Cloning 
Vector 连接 ,转化 到 Trans1-Tl 感受 态 细 胞 内 , 筛选 
阳性 单 克隆 进行 菌落 PCR 检测 ,之 后 送 到 上 海 生 工 
生物 工程 股份 有 限 公司 进行 双向 测序 。 
1.5 序列 整理 及 数据 分 析 

首先 将 正 反 向 引物 测 得 的 序列 导入 MEGA 7.0 
软件 (Kumar et al., 2016) 中 , 结合 峰 图 识别 软件 
Chromas 2. 6 校正 序列 ,再 进行 多 序列 比 对 。 利 用 
NCBI 数据 库 对 获得 的 序列 进行 BLAST 比 对 分 析 ， 
以 确定 扩 增 的 基因 序列 来 自 线粒体 COT 基因 。 

利用 MEGA 7.0 软件 (Kumar et al., 2016) X1 Z& 
PENE miCOI 基因 序列 进行 变异 分 析 , 计算 变异 位 
点 数 .简约 信息 位 点 数 、 碱 基 组 成 等 信息 。 同 时 , 基 
于 Kimura-2-Parameter( K2P ) 模 型 ,以 美洲 斑 潜 晶 工 . 
sativae PA X BIETER L. huidobrensis , — Wr TE EXE ES Wig 
L. trifolii, 38 Jp PEYS Wa L. bryoniae IUE 3 EZ TES Me 
Chromatomyia horticola 的 相应 序列 为 外 群 ,采用 非 
加 权 组 平均 (UPGMA ) 法 构建 单 倍 型 系统 发 育 树 , 系 
统 树 各 分 支 的 置信 和 度 (bootstrap ) 均 进行 1 000 KE 
复 检 验 。 利 用 DnaSP 6. 1 软件 ( Rozas et al., 2017) 
对 不 同 地 理 种 群 的 萤 斑 潜 蝇 进行 单 倍 型 多 样 性 
(Hd) AERE FEVE (Pi) fII Tajima’ s D 中 性 检验 等 
分 析 。 应 用 Network 4. 6 基于 Median-joining 算法 构 
建 单 倍 型 间 的 网 络 关系 图 。 












































应 用 Arlequin 3. 5 软件 (Excoffier and Lischer, 
2010) 计算 种 群 间 的 遗传 分 化 系数 (Fsr) (Hudson et 
al., 1992) .遗传 距离 [ Fs] (1 - Fer)]、 基 因 流 (Nm) 
等 参数 ;并 进行 分 子 方差 分 析 ( AMOVA) ,计算 遗传 
变异 在 种 群 间 和 种 群 内 的 方差 贡献 率 。 

根据 采集 地 经 纬度 计算 出 成 对 种 群 间 的 地 理 距 
离 , 同 时 结合 遗传 分 化 系数 (Fr) ,用 TFPGA 软件 
(Miller, 1997) rS Mantel 算法 检测 不 同 种 群 间 的 
遗传 距离 和 地 理 距离 的 相关 性 ,进行 10 000 次 置 
换 , 并 做 显著 性 检验 。 





























2 结果 


2.1 葱 斑 潜 蝇 线粒体 COI 基因 的 碱 基 组 成 及 单 倍 
型 分 布 

Tib B Z8 PESE P 12 个 地 理 种 群 共 253 条 
mtCOI 基因 序列 比 对 、 剪 切 后 得 到 759 bp 的 片段 ， 
分 析 表 明 ,无 碱 基 的 搬 和 人 或 缺失 , 具 明 显 的 A/T 碱 
基 偏 倚 性 , A + 了 平均 含量 为 70. 3% ;序列 共有 20 
个 变异 位 点 , 占 碱 基 总 数 的 2. 6496 ,包括 13 个 单 变 
异 多 态 性 位 点 和 7 个 简约 信息 位 点 。 

单 倍 型 的 数量 及 其 分 布 和 频次 的 结果 显示 ( 表 
2) , 供 试 种 群 中 共 发 现 13 种 单 倍 型 Hapl - Hap13 
( GenBank 登录 号 : MN700917 - MN700929) ,包括 7 
个 共享 单 倍 型 ( Hapl，Hap3，Hap6，Hap8 - 
Hap11 ) ,其 余 6 个 单 倍 型 为 相应 地 理 种 群 独 享 单 倍 
型 ,其 中 Hap2 属于 北京 海淀 种 群 (BJHD) ，Hap4 和 
Hap5 属于 甘肃 白银 种 群 (GSBY) ，Hap7 属于 海南 
陵 水 种 群 (HNLS ) Hapl2 和 Hap13 属于 吉林 通化 
种 群 (JLTH) 。Hapl 出 现 频率 最 高 , 占 总 样本 的 
81. 82% , 且 分 布 最 广 ,为 所 有 地 理 种 群 共有 ,其 余 
12 种 单 倍 型 出 现 频率 均 较 低 ( 和 3.95% ) 。 各 地 理 
种 群 中 ,吉林 通化 种 群 (JLTH) 的 单 倍 型 种 类 最 多 (5 
种 ) ,其 次 为 河北 石家庄 种 群 (HBSJZ) 河南 许昌 种 
群 (HNXC) 和 山东 济南 种 群 (SDJN) 有 4 种 ,而 河北 
保定 种 群 (HBBD) 与 内 蒙 乌 海 种 群 (NMWH) 仅 有 1 
种 单 倍 型 。 
2.2 葱 斑 潜 蝇 种 群 的 遗传 多 样 性 及 中 性 检验 

萄 斑 漆 蝇 各 种 群 mtCOI 基因 单 倍 型 多 样 性 
(Hd) 核 昔 酸 多 样 性 ( 户 )、 核 音 酸 平均 差异 数 ( 天 ) 、 
Tajima”s D 中 性 检测 分 析 结 果 显 示 ( 表 2) ,总 群体 
单 倍 型 多 样 性 ( Hd) 为 0.327, 核 芽 酸 多 样 性 ( Pi) 为 
0. 00159, 4X df R F Hg 23 96 XX (K) 为 
1.21011, 中 性 检验 结果 不 显著 (Tajima’sD = 
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A2 BE 12 个 不 同 地 理 种 群 基于 mtCOI 基因 序列 的 遗传 多 样 性 分 析 及 中 性 检验 
Table 2 Genetic diversity and neutrality test of Liriomyza chinensis among 12 geographical populations 


based on mtCOI gene sequence 























PERFIS 

单 倍 型 多 样 性 AERLE Tajima’ D 

单 倍 型 数 种 群 单 倍 型 组 成 (个 体 数 ) 差异 数 
种 群 li Haplotype Nucleotide did 中 性 检验 
Haplotype Haplotype composition within population Average number 
Population diversity diversity Tajima' D 
number ( number of individuals) of nucleotide 
Hd Pi neutrality test 
difference K i 
Hapl (8), Hap8 (2), Hapll (4), Hapl2 (5), 
JLTH 5 pot P 1 BA PV 0. 763 0.00716 5.43684 0. 50355 
Hap13 (1) 
BJHD 3 Hapl (19), Hap2 (1), Hap3 (1) 0.186 0.00132 1.00000 - 1.25759 
NMWH 1 Hapl (10) 0 0 - - 
HBBD 1 Hapl (20) 0 0 - - 

HBSJZ 4 Hapl (14), Hap8 (1), Hap9 (1) , HaplO (4) 0. 489 0.0015 1.13684 — 1.89783 * 
SDWF 2 Hapl (20), Hapl0 (1) 0.095 0. 00013 0.09524 -1.16356 
SDYT 3 Hapl (24), HaplO (1), Hapll (1) 0.151 0.00111 0. 84615 -2.34784 ** 
SDJN 4 Hapl (17), Hap3 (1), Hap6 (1) , HaplO (1) 0. 284 0.00092 0. 70000 -2.12144 ** 
SDLY 3 Hapl (16) , HaplO (3) , Hapll (1) 0.353 0.00167 1.26842 -2.09572** 
GSBY 3 Hapl (21), Hap4 (1), Hap5 (1) 0.170 0. 00023 0.17391 -1.51496 * 
HNXC 4 Hapl (15), Hap3 (1), Hap9 (6) , Hapll (2) 0.562 0.00315 2.38768 -1.43315* 
HNLS 3 Hapl (23), Hap6 (4) , Hap7 (2) 0.315 0. 00043 0. 32540 -0.76753 
总 计 Total 12 Hapl - Hap13 0. 327 0.00159 1.21011 - 1.64799 


Hapl - Hap13; 单 倍 型 Haplotypes. 单 倍 型 后 括号 内 的 数字 为 具有 该 单 倍 型 的 个 体 数 。The numbers in brackets after haplotype are the number of 


individuals with this haplotype. * P «0.05; ** P <0.01. 








-1.64799, P 20.05) ,说 明 供 试 萄 斑 潜 蝇 总 群体 在 ”与 其 
过 去 较 近 的 一 段 历 史 时 期 内 未 经 历 明显 的 群体 扩张 





二 他 单 倍 型 遗传 距离 最 大 ,遗传 距离 达到 了 
0. 012 ~0. 020; 其 次 是 单 倍 型 Hap2, Hap3 Hap8 和 


























和 持续 增长 模式 ,群体 大 小 保持 相对 稳定 状态 。 各 
地 理 种 群 单 倍 型 多 样 性 (Hd) 在 0 ~0.763 [HR], zd 
酸 多 样 性 (Pi) 在 0~0.00716 H , 核 苷 酸 平均 差异 数 
在 0. 09524 ~ 5. 43684 间 ,吉林 通化 种 群 (JLTH ) 的 
单 倍 型 多 样 性 最 高 (Hd 20.763) ,其 次 是 河南 许昌 
种 群 ( HNXC，Hd = 0. 562) 和 河北 石家庄 种 群 
(HBSJZ, Hd = 0. 489) ,而 河北 保定 种 群 ( HBBD, 
Hd =0) 和 内 蒙 乌 海 种 群 ( NMWH, Hd 20) Rf Ze FE 
性 最 低 。 各 种 群 间 的 Tajima’ s D 中 性 检测 分 析 结 
果 表 明 , 甘 肃 白银 种 群 (GSBY)、 河 北 石家庄 种 群 
( HBSJZ) 和 河南 许昌 种 群 (HNXC ) 表现 显著 ,山东 










































































Hap13 ,这 4 个 单 倍 型 与 其 他 单 倍 型 之 间 的 遗传 距 
离 在 0.003 ~0.019 之 间 ; 其 余 7 个 单 倍 型 之 间 遗 传 
距离 均 很 接近 ( <0.003)。 

从 萄 斑 潜 蝇 mtC01 基因 单 倍 型 系统 发 生 树 可 
见 ( 图 1) , ZR DIETE NS 13 种 单 倍 型 中 , Hapll 单独 形 
成 一 文系 ,其 他 12 种 单 倍 型 聚 为 另 一 文系 。 这 12 
种 单 倍 型 聚 成 的 大 支 又 可 分 为 两 小 支 系 ,第 一 小 文 
包含 Hap2, Hap3, Hap8 和 Hap13 共 4 个 单 倍 型 ,第 
二 小 支 系 包含 Hapl, Hap4 - Hap7, Hap9 - 10 及 
Hapl2 dt 8 种 单 倍 型 。 

为 了 进一步 直观 体现 出 萄 斑 潜 蝇 mtCOI 基因 









































济南 种 群 (SDJN) .山东 临沂 种 群 





(SDLY ) 和 山东 烟 


单 们 型 之 间 的 遗传 分 化 及 在 种 群 





中 的 分 布 频率 , 构 





台 种 群 (SDYT) 表现 极 显著 ,说明 这 6 T DC ZB 
洲 晶 种 群 在 较 近 的 历史 时 期 经 历 了 种 群 扩 张 , 而 其 
他 地 区 的 群体 表现 均 不 显著 。 
2.3 PEBH mtCOI 基因 单 倍 型 遗 
进化 分 析 

遗传 分 化 的 结果 显示 ( 表 3) ,各 单 倍 型 间 的 
K2P 遗传 距离 范围 为 0.001 ~ 0. 020 , 单 倍 型 Hapll 






































传 距离 及 系统 























建 13 种 单 倍 型 的 Network 网 络 关 系 图 (图 2) 。 结 果 
显示 ,以 单 倍 型 Hapl 为 中 心 ,其 他 单 倍 型 均 围 绕 这 
一 单 倍 型 呈 星 形 辐射 分 布 状 , 据 此 推测 分 布 广泛 且 
分 布 频次 最 高 的 Hapl 为 祖先 单 倍 型 。 
2.4 不 同 地 理 种 群 间 的 遗传 分 化 .基因 交流 和 分 子 
方差 分 析 

萄 斑 潜 蝇 总 群体 的 基因 流 (Nm) 为 3.33629, 遗 
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表 3 基于 葱 斑 潜 蝇 地 理 种 群 mtCOI 基因 序列 不 同 单 倍 型 的 遗传 距离 


Table 3 Genetic distances of different haplotypes of Liriomyza chinensis geographical populations 


based on mtCOI gene sequence 


单 倍 型 Haplotype ^ Hapl Hap2 Hap3 Hap4 Hap5 Hap6 Hap7 Hap8 


Hap9 Hapl0  Hapll Hapl2  Hapl3 





















































Hapl 0.003 0.003  Á 0.001 40.001 40.001 0.001 0.003 00.001 0.001 0.004 0.001 0.003 
Hap2 0.008 0.001 0.003 0.003 0.003 0.003 0.002 à 0.003 0.004 0.005 0.003  Á 0.002 
Hap3 0.007 0.001 0.003 0.003 0.003 0.003 0.002 à 0.003 0.003 0.005 0.003 0.001 
Hap4 0.001 0.007 0.005 0.002 0.002 0.002 0.003 0.002 à 0.002 à 0.004  Á 0.002 0.003 
Hap5 0.001 0.007 0.008  Á 0.003 0.002 0.002 0.004 0.002 0.002 4 0.004 0.002 0.003 
Hap6 0.001 0.009 0.008 0.003  Á 0.003 0.002 0.004 0.002 0.002 0.005  Á 0.002 0.003 
Hap7 0.001 0.007 0.005 0.003 0.003 0.003 0.003 0.002 0.002 0.005 0.002 0.003 
Hap8 0.009 0.004 0.003 0.008 0.011 0.011 0.008 0.004 0.004  Á 0.005 0.004 0.001 
Hap9 0.001 0.009 0.008 0.003 0.003 0.003 0.003 0.011 0.002 0.004 0.002 0.003 
Hap10 0.001 0.009 0.008 0.003 0.003 0.003 0.003 0.011 0.003 0.005 0.002 0.004 
Hapll 0.013 0.016 0.017 0.012 0.012 0.015 0.015 0.020 0.015 0.015 0.004 0.005 
Hap12 0.001 0.009 0.008 0.003 0.003 0.003 0.003 0.011 0.003 0.003 0.015 0. 004 
Hap13 0.008 0.003 0.001 0.007 0.009 0.009 0.007 0.001 0.009 0.009 0.019 0.009 
对 角 线 以 下 为 遗传 距离 平均 值 ,对 角 线 以 上 为 标准 误 。Data below the diagonal are means of the genetic distance and those above the diagonal are the 
standard error. 
[ Hapl 
Hap9 
| Hap7 
Hap6 
| Hap12 
76 
Hap10 
| Hap4 
87 — Hap5 
Hap2 
100 Hap3 
82 | Hap8 
Hapl13 
Hapll 
VRE Liriomyza sativae (AB721353.1) 
98 O EE Liriomyza trifolii (AB721357.1) 
番茄 斑 潜 蝎 Liriomyza bryoniae (AB721335.1) 
99 RRES Liriomyza huidobrensis (AB721351.1) 
Wü c EB Chromatomyia horticola (AB721328.1) 
te — — —— E ER oa E 
0.080 0.060 0.040 0.020 0.000 





图 1 UPGMA 法 构建 的 基于 mCOI 基因 序列 的 葱 斑 潜 蝇 地 理 种 群 


及 近 缘 种 的 单 倍 型 系统 发 育 树 








Fig. 1 Phylogenetic tree of haplotypes of Liriomyza chinensis geographical populations and 


related species constructed by UPGMA method based on mtCOI 
Hapl - Hap13; 单 倍 型 Haplotypes. 





gene sequence 
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Ho 
图 2 JEDER TE mtCOI 基因 单 倍 型 中 介 网 络 关系 图 
Fig. 2 Median-joining haplotype network of Liriomyza chinensis geographical populations based on mtCOI gene sequence 

每 个 圆 形 表示 一 种 单 倍 型 , 圆 形 面积 反映 出 其 单 倍 型 在 总 群体 中 的 出 现 频率 ,H_1 - H_13 分 别 与 单 倍 型 Hapl - Hap13 相对 应 ,mvl 和 m2 表 














示 可 能 的 突变 位 点 。Each circle indicates one haplotype, and the circle size represents the frequency of haplotypes appearing in all populations tested. 


H 1 - H. 13 correspond to haplotype Hapl - Hap13, respectively, and mvl and mv2 indicate possible mutation sites. 


传 分 化 系数 (Fsr) 为 0. 06971, 符合 1 和 Nm «4, 
0.05 < For <0. 15 ,说 明 群 体 间 存在 中 等 程度 的 遗传 
分 化 ,基因 交流 较 充 分 ( Wright，1951 ) 。 各 种 群 间 
的 Fsr 在 -0.045 ~0.191 之 间 ( 表 4)。 吉 林 通 化 种 
群 (JLTH) 与 甘肃 白银 种 群 (CSBY) ` 海 南陵 水 种 群 
(HNIS) 河北 保定 种 群 (HBBD)、 山 东 潍 坊 种 秋 
(SDWF) 间 的 遗传 分 化 程度 较 大 (0. 15 < For < 
0. 25) ,同时 基因 交流 较为 充分 (1 和 Mn <4)( 表 4) 。 

XAET Z8. DETEE NB Fb EE fT 09 ^r T 27 25 2T Br 
(AMOVA ) HR CX 5 ) RHH , ADEE BAPE N Pe E 
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TH 




















异 显著 (已 <0.05 ) , 占 总 变异 的 93. 23% , 说明 遗 传 
变异 主要 来 自 于 种 群 内 部 ;种 群 间 遗传 变异 水 平 低 ， 
仅 占 总 变异 的 6. 97%; 不 同 寄 主 间 (Fo = 
- 0.00485 <0. 05 ) 表现 出 较 低 的 遗传 分 化 水 平 ,各 
寄主 植物 间 的 遗传 变异 不 明显 (已 >0.05) ;同一 寄 
主 不 同 地 理 种 群 间 (Fsc =0.07216 ) 表 现 出 中 等 程度 
的 遗传 分 化 ,遗传 变异 显著 (P<0.05)。 

对 所 有 种 群 进行 Mantel 检验 ,发现 遗传 距离 与 
地 理 距离 之 间 不 存在 显著 的 相关 性 (r = 0. 3820, 
已 >0.05)( 图 3)。 
























































ARA 藻 斑 潜 蝇 成 对 种 群 间 遗 传 分 化 系数 (Fsr, 下 三 角 ) 与 种 群 间 的 基因 流 (Nm, 上 三 角 ) 


Table 4 Pairwise genetic differentiation coefficient (Fs, below the diagonal) and gene flow ( Nm, above the diagonal) 


between geographical populations of Liriomyza chinensis 





nm JLTH BJHD NMWH HBBD HBSJZ SDWF SDYT SDJN SDLY GSBY HNXC HNLS 

JLTH 4.381 4. 720 2.557 3.708 2.507 3.596 3.514 5.044 2.606 9.192 2.122 

BJHD 0.102 inf 11.351 67.200 11.500 33.770 inf 22.050 21.589 16.157 8.281 

NMWH 0.096 -0.005 inf 180. 000 inf inf inf inf inf 32.991 23.817 
HBBD 0. 164 0. 042 0. 000 9. 000 inf inf inf 20.083 inf 7.460 7.075 

HBSJZ 0.119 0.007 0.003 0.053 21.832 34.547 inf inf 10.184 13.031 6.632 

SDWF 0.166 0.042 -0.041 -0.002 0.022 inf inf 113.546 inf 7.2773 7.858 

SDYT 0.122 0.015 -0.045 -0.010 0.014 -0.014 inf inf inf 23.983 17.340 
SDJN 0.125 -0.029 -0.040 0.000 -0.015 -0.010 -0.011 inf inf 13.214 320.178 
SDLY 0.090 0.022 -0.020 0.024 -0.015 0.004 -0.032  -0.004 31.680 44.078 9.2774 

GSBY 0.161 0.023 -0.043 -0.006 0.047 -0.001 -0.016 -0.006 0.016 8.405 8.426 

HNXC 0.052 0.030 0.015 0.063 0.037 0.064 0.020 0.036 0.011 0.056 5.445 

HNLS 0.191 0.057 0.021 0.066 0.070 0. 060 0.028 0.002 0.051 0.056 0.084 


inf: 无 限 大 Infinite. 
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Table 5 Analysis of molecular variance ( AMOVA) within and among populations 
of Liriomyza chinensis based on mtCOI gene 
; " 平方 和 方差 组 成 变异 百分比 D 

变异 来 源 自由 度 . 固定 系数 

Sum of Variance Percentage of 
Source of variation d. f. Fixation index 

squares components variation 
种 群 间 Among populations 11 16.030 0.042Va 6.97 
寄主 植物 间 Among host plant species 1 1.156 — 0. 00295 Va -0.49 Fer = -0.00485, P »0.05 
寄主 内 地 理 种 群 间 

10 14. 874 0.04403 Vb 7.25 F4; 20.07216, P «0.05 

Among geographic populations within host plant species 
种 群 内 Within population 241 136. 444 0. 56616 Vc 93.23 Fsr =0.06766, P «0.05 











Va, Vb, Ve: 分 别 表示 方差 组 分 的 数量 Number of variance components; For: 寄主 种 群 间 变 异 Variation among host populations; Fsc: 寄主 种 群 内 


变异 Variation within host populations; Fop: 种 群 内 变异 Variation within a population. 
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图 3 DEVENE 12 个 地 理 种 群 间 遗 传 距离 与 地 理 距 离 的 相关 性 


Fig. 3 Correlation between the genetic distance and the geographical distance 














for 12 geographical populations of Liriomyza chinensis 


Fgr/ (Y -了 sz) 表示 遗传 距离 。The genetic distance is expressed as Fsr/(1 — Fs). 


线粒体 DNA ( mtDNA ) 是 常用 的 分 子 遗 传 标记 ， 
其 中 的 线粒体 细胞 色素 氧化 酶 亚 基 I(COI) 基因 在 
昆虫 种 类 鉴定 、 隐 存 种 的 发 现 、 种 群 遗传 多 样 性 MJ 
种 的 系统 进化 等 研究 方面 均 得 到 广泛 应 用 ( DeSalle 
et al., 1987 ; 成 新 跃 等 , 2000; Hebert et al., 2003; 
何 恒 果 ，2008; Norgate et al., 2009; Ye et al., 
2018), ÆF mtCOI 基因 ,本 研究 对 我 国 8 个 省 / 直 
辖 市 的 12 个 不 同 地 区 的 和 瓯 斑 潜 蝇 种 群 的 遗传 多 样 
性 、 遗 传 分 化 程度 基因 交流 等 进行 了 深入 分 析 。 

本 研究 结果 表明 , 莹 斑 潜 蝇 mtDNA. COL 基因 序 








Wang et al., 2014; Chen et al., 2018) 相符 。 在 12 
个 地 理 种 群 中 共 定 义 了 13 个 单 倍 型 ,包括 7 个 共享 
单 倍 型 和 6 个 独 享 单 倍 型 ,说 明 不 同 地 理 区 域 的 葱 
斑 潜 晶 在 保持 着 稳定 遗传 变异 的 同时 ,也 产生 了 适 
应 特定 地 区 环境 的 遗传 变异 。 在 7 个 共享 单 倍 型 
中 , 单 倍 型 Hapl 在 所 有 地 理 种 群 中 均 有 分 布 且 出 
现 频率 最 高 ,同时 结合 单 倍 型 网 络 结构 图 来 看 ,可 认 
为 单 倍 型 Hapl 为 葱 斑 潜 晶 群 体 中 适应 性 较 强 并 广 
泛 分 布 的 原始 单 倍 型 。 

本 研究 中 欧 斑 潜 临 总 群体 的 单 倍 型 多 样 性 较 低 
( 表 2)。 有 人 研究 报道 ,地 理 距离 .寄主 植物 种 类 以 及 
生境 等 是 影响 昆虫 的 种 群 遗 传 多 样 性 和 遗传 结构 的 
重要 因素 (Morgan et al., 2000; Scheffer and Lewis, 




















列 片段 的 碱 基 组 成 表现 较 强 的 A/T 偏 倚 性 ,与 典型 
的 昆虫 线粒体 碱 基 组 成 特点 (Simon et al., 1994; 


2006; Tuda et al., 2014 )。 而 本 研究 Mantel 检验 结 
果 显 示 , 地 理 距离 与 欧 斑 潜 蝇 的 遗传 距离 间 没 有 相 
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关 性 , 即 地 理 距 离 并 不 影响 葱 斑 潜 蝇 的 遗传 分 化 
(图 3)。 从 寄主 植物 的 角度 来 讲 , 敬 斑 洲 晶 是 寡 食 
性 昆虫 , 仅 取 食 百 合 科 萄 属 , 旦 主 要 寄主 植物 为 大 
萄 .小 萄 ,在 其 他 植物 上 少见 发 生 , 种 类 单一 的 寄主 
植物 可 能 在 一 定 程度 上 限制 了 葱 斑 潜 蝇 种 群 的 遗传 
分 化 。 在 12 个 地 理 种 群 中 , 仪 吉 林 通 化 种 群 的 单 倍 
型 多 样 性 表现 较 高 (Hd =0.763)( 表 2) 。 结 合 采样 
地 特征 分 析 其 原因 ,吉林 通化 的 萄 斑 潜 蝇 为 农家 自 
留 地 采集 ,与 商业 化 大 规模 种 植 的 萄 相 比 ,化 学 农药 
的 使 用 量 较 少 , 较 低 的 药剂 选择 压力 可 能 增加 了 区 
斑 潜 蝇 的 遗传 变异 。 总 的 来 说 ,作者 认为 萄 斑 潜 晶 
在 吉林 通化 可 能 受到 人 工 选择 压力 较 小 以 及 其 他 原 
因 ,造成 了 其 较 高 的 遗传 多 样 性 。 

遗传 分 化 系数 (Fsr) 和 基因 流 (Nm) 是 研究 种 
群 遗传 分 化 的 重要 参数 ,可 以 解释 影响 种 群发 生 遗 
传 分 化 的 因素 (Rousset，1997) 。Tajima' s D 检验 可 
揭示 物种 在 一 定 历史 时 期 内 的 种 群 动态 等 情况 , 通 
过 检验 值 来 判断 种 群 是 否 存在 扩张 及 瓶颈 效应 
(Tajima, 1989) 。 在 本 研究 中 , A DEWE IE A REA KY 
遗传 分 化 系数 (Fir) 为 0. 06971 (0.05 < For < 
0. 15) ,遗传 分 化 程度 达到 中 等 水 平 ,遗传 变异 主要 
来 自 于 种 群 内 ; 基因 流水 平 (Nm) 为 3. 33629 (1 < 
Nm <4) ,群体 交流 较为 充分 ( 表 4) ,说 明和 葱 斑 潜 晶 
有 某 种 渠道 可 以 进行 基因 交流 。 基 因 交 流水 平 可 能 
与 昆虫 自身 的 迁徙 能 力 交通 运输 气流 扩散 等 因素 
密切 相关 。 而 斑 潜 蝇 属 害 忠 普遍 飞行 能 力 较 弱 , 据 
报道 ,美洲 斑 潜 蝇 平均 飞行 距离 仅 为 0.95 km( 雷 仲 
仁 等 , 2002) 。 萄 斑 潜 晶 与 美洲 斑 潜 蝇 体 型 相近 ,为 
韭 迁 飞 性 害虫 ,故而 推测 萤 斑 潜 蝇 通过 飞行 进行 不 
同 地 区 基因 交流 的 可 能 性 较 小 。 随 着 物流 行业 的 发 
展 ,交通 运输 也 逐渐 成 为 昆虫 进行 长 距离 基因 交流 
的 一 种 重要 途径 (Tuda et al., 2014). Z DEVS WE fE 
害 的 主要 寄主 植物 是 大 萄 ,而 我 国 的 大 莹 种 植 主要 
集中 在 山东 河北、 河南 所 在 的 第 一 大 主 产 区 ,频繁 
的 交通 运输 可 能 是 敬 斑 潜 晶 群体 基因 交流 较 充 分 的 
重要 原因 ,而 葱 斑 潜 蝇 是 否 存在 气流 扩散 方式 ,这 方 
面 的 研究 还 未 有 报道 。 中 性 检验 结果 表明 , 萄 斑 潜 
晶 总 群体 在 较 近 的 历史 时 期 内 未 经 历 明 显 的 种 群 扩 
张 或 种 群 增长 ,总 群体 保持 稳定 状态 。 

本 研究 仅 以 mtDNA. COI 基因 长 为 759 bp 的 单 
个 片段 作为 分 子 遗 传 标记 ,为 研究 葱 斑 潜 蝇 的 系统 
进化 等 提供 了 基础 性 遗传 资料 。 本 研究 发 现 新 斑 洪 
晶 的 遗传 多 样 性 水 平 低 , 后 续 还 需要 增加 不 同 生境 
类 型 .寄主 植物 种 类 、 虫 源 采集 量 等 数据 ,才能 更 深 
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